The isoscalar monopole resonance of the alpha particle: 
a prism to nuclear Hamiltonians 
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We present an ab-initio study of the isoscalar monopole excitations of 4 He using different realistic 
nuclear interactions, including modern effective field theory potentials. In particular we concentrate 
on the transition form factor Fm to the narrow + resonance close to threshold. Fm exhibits a 
strong potential model dependence, and can serve as a kind of prism to distinguish among different 
nuclear force models. Comparing to the measurements obtained from inelastic electron scattering 
off 4 Hc, one finds that the state-of-the-art theoretical transition form factors are at variance with 
experimental data, especially in the case of effective field theory potentials. We discuss some possible 
reasons for such discrepancy, which still remains a puzzle. 
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The isoscalar monopole strength of large nuclei has 
been extensively studied since the discovery of a giant 
monopole resonance in 144 Sm and 208 Pb [lj. The rea- 
son for the great interest in such excitations originates 
from their connection to the incompressibility modulus 
of infinite nuclear matter 0, H| . The alpha particle is a 
light nucleus, that however has a binding energy per par- 
ticle similar to that of large systems and a high central 
density. While it possesses no bound excited states, it ex- 
hibits a very pronounced narrow resonance ( 4 He*) with 
the same quantum numbers + as the ground state, i.e., 
an isoscalar monopole resonance. Today, the develop- 
ment of few-body theories has reached a point, where an 
ab-initio calculation of the four-body isoscalar monopole 
transition strength can be carried out with high preci- 
sion. As will become evident in the following, the com- 
parison of such four-body results with experimental data 
can serve as a stringent test for nuclear Hamiltonians, 
that are the sole ingredients of an ab-initio quantum me- 
chanical approach. 

The four-nuclcon isoscalar monopole resonance is lo- 
cated at E e * p = -8.20 ± 0.05 MeV, with a width of 
270±50 keV just above the first two-body break-up 
threshold Ef hr = —8.48 MeV into a proton and a triton 
and below the next threshold E\) lr = -7.74 MeV into a 
neutron and a 3 Hc. A summary about the experimental 
studies of the spectrum of 4 He can be found in Ref. Q. 
Valuable information about the nature of the resonance 
is given by the transition form factor Fm (<z) measured in 
electron scattering experiments ( 4 He(e, e') 4 He*) at vari- 
ous momentum transfer q. Similarly to the case of the 
elastic form factor, the q dependence of Fm reflects the 
dynamics at various interaction ranges. 

The progress in ab-initio few-body methods allows to- 
day to obtain accurate results for observables in light 
nuclear systems using realistic potential models (see re- 
view Q). In recent years the debate regarding poten- 
tial models has boosted, especially after the introduction 
of the effective field theory (EFT) strategy in nuclear 
physics Q- At present, both phenomenological realistic 



and chiral EFT potentials are used in ab-initio calcula- 
tions, but only for very few observables large differences 
are found, e.g., for the polarization observable A y of p- 
3 He scattering Q. In this letter, we point out that the 
calculated 4 He isoscalar monopole resonance transition 
form factor Fm(q) depends dramatically on the nuclear 
Hamiltonian. Thus, it can serve as a kind of prism to 
distinguish among nuclear force models. 

Main Results. The isoscalar monopole strength 
Sm (<7j til) is in general a function of q and the energy 
transfer uj. It is given by 



S M {q,u) = dn\(n\M(q)\0)\ 2 5(uj - E n + E ) (I) 

-M(q)\0), 



■ l/ ro <0|Ait (g) _^_ 



where |0),|n) and E ,E n are eigenfunctions and eigen- 
values of the nuclear Hamiltonian H , and 



M(q) = 



G s E (q) 



(2) 



is the isoscalar monopole operator. Here G s E (q) — 
G E (a) + G E (q) is the nucleon electric isoscalar form fac- 
tor [9(, Ti is the nuclcon's position, and jo is the spherical 
Bessel function of th order. The monopole strength can 
be written as a sum of a resonance term and a non- 
resonant background contribution S^, 



S M (q,u) = S r M(q,u) + S^(q,uj) 



(3) 



For a narrow resonance one defines the resonance transi- 
tion form factor 



\F M (q)\ 2 = ± 



(q,u) 



(4) 



In Fig. [TJ we show results for Fm (<?) with two differ- 
ent Hamiltonians including realistic three-nucleon forces 
(3NF) in comparison to experimental data from inelas- 
tic electron scattering 
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1. (Color online) Theoretical transition form factor 
q' 2 )\ 2 with G% — calculated with various force mod- 
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FIG. 2. (Color online) Elastic form factor \F e ,(q 2 )\ of 4 He 
calculated with various force models: AV18+UIX (full line); 



els: AV18+UIX (full line), N 3 LO+N 2 LO (red band); re- N 3 LO+N 2 LO (red band); result fromjlj with AV18+UIX 



suit from [l^l (dot-dashed). Data from Frosch et al. [1 
Walcher Hi and Kobschall et al. GJ. 



(dot-dashed). Data from Frosch et al. 



use (i) the AV18 [13| NN potential plus the UIX (14 1 
3NF, (ii) an EFT based potential, where we take the NN 
potential [l5[ at fourth order (N 3 LO) in the chiral ex- 
pansion augmented by a 3NF at order N 2 LO [l6[. The 
Coulomb potential is taken into account in all calcula- 
tions. Both the EFT and the AV18 NN potentials repro- 
duce the NN scattering phase shifts with high precision 
(% 2 /datum rj 1). In the EFT calculations, two different 
paramcterizations of the 3NF have been used, leading to 
the red band in Fig. [T] The chiral low energy constants 
Cd and ce have been determined either by setting Co 
to a reasonable value and then fitting ce to the three- 
nucleon binding energies [l6[ (cd = 1 and ce = —0.029) 
or by fitting to the 3 H binding energy and beta decay [l7( 
(cd = —0.2 and ce = —0.205). We also display the re- 
sult of a previous calculation by Hiyama et al. [lfj , with 
the AV8' potential, a reduced version of AV18, and a 
simplified central 3NF, fitted to the binding energy of 
3 H. All three Hamiltonians reproduce the 4 He experimen- 
tal binding energy within one percent. Surprisingly, the 
results for Fj^[{q) strongly depend on the Hamiltonian. 
Furthermore, the realistic Hamiltonians fail to reproduce 
the experimental data. In particular, this is true for the 
EFT forces that predict a transition form factor twice as 
large as the measured one. 

In contrast, the realistic Hamiltonians lead to rather 
similar results for the elastic form factor F e i(q) of 4 He, 
defined as 



F el (q) = -(0\M(q)\0) 



(5) 



In Fig. [2j F e i(q) is shown for the AV18+UIX model and 
for the chiral EFT potentials. The fact that the results 
do not differ significantly is not very surprising, since 
both Hamiltonians give a very similar result for the ra- 



TABLE I. Ground state energies in MeV for 3 H, 3 He and 
4 He with N 3 LO [I| and N 2 LO (parameterizations from []]| 
or [I3|). Comparison of present results (EIHH) with no core 
shell model (NCSM) and hyperspherical harmonics (HH) re- 
sults. 



3NF from [16] 


EIHH 


NCSM [16] HH [25] 


Nature 


3 H 


-8.474(1) 


-8.473(5) -8.474 


-8.48 


3 He 


-7.734(1) 


-7.733 


-7.72 


4 He 


-28.357(7) 


-28.34(2) -28.37 


-28.30 


3NF from [17] 


EIHH 


NCSM [17] 


Nature 


3 H 


-8.472(3) 


-8.473(4) 


-8.48 


3 He 


-7.727(4) 


-7.727(4) 


-7.72 


4 He 


-28.507(7) 


-28.50(2) 


-28.30 



dius: 1.432(2) fm [20| for AV18+UIX and 1.464(2) fm 
for N 3 LO plus the N 2 LO of [T^ which is not far from 
the experimental value of 1.463(6) fm (obtained from 
the charge radius of Ref. [2l| as explained in (2p|). Also 
shown in Fig. [5] is the result by Viviani et al. [18[ with 
the AV18+UIX potential, which is indistinguishable from 
ours, proving the level of accuracy of contemporary four- 
body calculations. 

Calculational Method. Our calculations are carried 
out using the effective interaction hyperspherical har- 
monics method (EIHH) 0, H3]. In Tabic [fl wc present 
the EIHH binding energies of three- and four-body nuclei 
obtained from EFT potentials including 3NF. Compar- 
ing our results to other calculations we find agreement at 
the 10 keV level. Results for Sm (q> w ) are often obtained 
by discrctizing the continuum. In this case the Hamilto- 
nian is represented on a finite basis of square integrablc 
functions and is then diagonalized to obtain eigenvalues 
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e v and eigcnfunctions \v). In this way one achieves an 
ill defined discretized representation of S.m(q,w)- In our 
Lorentz integral transform (LIT) approach [2(| [2?| the 
continuum discretization is used to reach the continuum 
spectrum as follows. We calculate 



C M (q,a,T) = —Im(0\M^q) 



1 



a -H + Eo + iT 



M(q)\0), 
(6) 

where T is finite. It is easy to proove that Cm (q, a, T) is 
connected to Sm (<Z, w) by an integral transform with a 
Lorentzian kernel K(uj,a, T) 



ir (u}+E -cr) 2 +r 2 ' 



C M (q, <r, T) = J K(w, a, T) S M 



(7) 



Since T is finite Cm (<?, o - , F) can be calculated with 
bound-state methods, namely it is legitimate to represent 
the Hamiltonian on a basis of square integrable functions. 
After having reached convergence in the energy region of 
interest, with N basis functions, the LIT of SVi(<Z,w) is 
given by 



N 



C M (q,<7,T) = -V 

7T ' » 



IH^(9)|0>| 



=J (a-e,+.Bo) 2 + r 2 



(8) 



In the next step an inversion of the transform ([7]) has to 
be performed to determine Sm(q,w)- This same proce- 
dure has been followed in Ref. [28|, |29[ but for the com- 
plete nuclear charge operator. There, however, the reso- 
nance region was not considered. 

The diagonalization of the Hamiltonian allows us 
to identify the resonance energy at Er = —7.40(20) 
McV for AV18+UIX and E R = -7.50(30) McV for 
N 3 LO+N 2 LO (extrapolated value). In fact wc find that 
for e v = E R the strength Sr = \(vn\M(q)\0)\ 2 is very 
pronounced. The strength Sr is composed of a transi- 
tion to the resonance and a background contribution. As 
explained in the following, by a proper inversion of our 
convergent results for £m(<1, cr, L) at T = 5 MeV we have 
been able to separate both contributions, i.e. FM(q) and 
Sm^u). 

For the inversion of the LIT we have used a regu- 
larization procedure (30| . which usually consists in ex- 
panding the response on a set of continuous functions 
and in fitting the calculated Cm (<3S f i T) on the cor- 
responding linear combination of the transformed ba- 
sis functions. In the presence of a narrow resonance 
an exp licit resonance form must be added to the ba- 
3l|, e.g., a Lorentzian with width 7 located at Er, 

; ) = 2^ (^E^M^hm? ■ In the ideal case both 

the resonance strength and its width can be resolved in 
the inversion process. However, in our present case, one 
finds that the best fits are obtained with a vanishing 
width, reflecting the absence of states \ v) in the vicinity of 
Er. Therefore, we treat the resonance as a S— function, 
which wc simulate numerically by setting 7 to a very 
small value. 
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FIG. 3. (Color online) (a) <S^(g,w) at q = 1.5 far 1 for 
AV18+UIX obtained with different values of fa (see text): 
f R = \FM{q)\ 2 (fall line), f R = 0.99|FM(g)| 2 (dashed 
line); fn = 1.01\Fm (?)| 2 (dot-dashed line); (b) non-resonant 
background contribution Sj^(q,bj): AV18+UIX (fall line); 
N 2 LO+N 3 LO (red band). 



Since in our calculation the resonance turns out to be 
a 5— function, the background LIT has the following form 



Am^o-,!") = -Cm (<?, er, F) 



fn(q) 



ir(a-E R + E )2+r 2 



(9) 

where f R remains to be determined in the inversion. 
Varying Jr one finds inversion results for S M that con- 
tain different resonance contributions. Results for vari- 
ous values of fn arc shown in Fig. [3^,. The case where 
the resonance contribution vanishes corresponds to the 
correct value of the transition form factor \FM(q)\ 2 = 
fR(q)/Z 2 and the obtained strength is just Sj^ (q, u) . Wc 
would like to emphasize that the results obtained are al- 
most independent of the assumed resonance width 7, so 
long it remains sufficiently small. For the AV18+UIX 
potential the relative size of the background reduc- 
tion, about 8%, is roughly q independent. For the 
N 3 LO+N 2 LO interactions the reduction varies between 
13% for q = 0.25 fnT 1 and 22% for q = 2 fin -1 . 

In Fig. [3}d, the non-resonant monopole strength S m is 
shown on a larger energy range, into the far continuum 
region. One sees quite a difference between the results 
with the EFT forces and the AV18+UIX potential. The 
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TABLE II. |F ol | and S R = \{u R \M(q)\0)\ 2 for q = 1.01 fin -1 
as a function of the grandangular momentum .ff max with 
N 3 LO+N 2 LO [H. 



12 



14 



16 



18 



\F ci 



0.6248 0.6244 0.6242 0.6241 



W 4 S r /4ttZ 2 4.59 4.75 4.85 4.87 



former lead to a lower low-energy peak and tail than the 
latter. These results show the power of the LIT approach, 
which enables one to calclate the strength in the far four- 
body continuum by reducing a scattering-state problem 
to a bound-state problem in a rigorous way. 

Analysis of the Results. The main findings of this Let- 
ter are the dramatic sensitivity of Fj^ (q) to the nuclear 
Hamiltonian and the large deviations of realistic calcula- 
tions from the available experimental data. These are un- 
expected results, so in the following we will try to analyze 
the possible sources for this puzzling situation. In Fig. Q] 
we presented all the available experimental data from 
Frosch et al. [TO], Walcher [|] and Kobschall et al. [13]. 
Fm{q) has been obtained in inelastic electron scatter- 
ing experiments measuring the longitudinal cross section 
at fixed values of q in the resonance region, subtracting 
the background contributions and integrating over the 
resonance structure. Even though one can contemplate 
the possibility of systematic experimental errors, the fact 
that the results correspond to three different sets of data, 
makes it less likely, so we will now discuss possible theo- 
retical explanations. 

(i) Our calculations are well converged in terms of the 
expansion of the wave functions. To elucidate this nu- 
merical aspect we show in Table |TT] both F c \ and Sr at 
q = 1.01 fm _1 for different grandangular momentum val- 
ues K„ 



which is a measure of the wave function ex- 
pansion in the EIHH method. Our results for Sr (and 
thus for -Fm) are converged at a few per cent level for all 
the implemented potentials, thus large numerical inaccu- 
racies can be ruled out. 

(ii) The isoscalar monopole operator in Eq. ^ is ob- 
tained from a multipole expansion of the usual one-body 
charge operator. The neglect of higher order operators 
(two- and three-body) is known as impulse approxima- 
tion. The impulse approximation describes the elastic 
form factor very well up to a momentum transfer of q = 2 
fur 1 (see Fig. For the AV18+UIX potential it has 



been shown explicitly by Viviani et al. 



that 



two-body charge operators are negligible to F e \(q) below 
q = 2 fin . This is not surprising since it is well known 
that two-body terms in the nuclear charge operator are 
of higher relativistic order. In fact, in EFT they appear 
only at fourth order [32[ • Thus, they are unlikely to be 
the source of the large discrepancy between theory and 
experiment. 



(Hi) The EFT potentials we have used are a NN inter- 
action at N 3 LO and 3NFs at N 2 LO. One may wonder 
whether N 3 LO corrections to the 3NF can have a large 
effect. However, we notice that the overall effect of the 
3NF on F M {q) is rather mild (about 10%). 
(iv) In our calculation the values of Er given by both 
potential models (AV18+UIX, N 3 LO+N 2 LO) lie beyond 
the two thresholds and about 700 keV above the exper- 
imental result. On the other hand, the simplified force 
model used by Hiyama et al. [l(| reproduces the experi- 
mental resonance position to within 100 keV, and at the 
same time predicts a smaller Fm(q), closer to the exper- 
imental data. One can envisage a correlation between 
the ability of a model to reproduce Er and Fm- In fact, 
if one considers that the transition form factor is the 
Fourier transform of the transition density from 4 Hc to 
He*, one can imagine that small differences in Er are 
reflected in the resonant wave functions and yield larger 
differences in the transition density. Similar conclusions 
have been drawn in Ref. [33| in the study of the dif- 
ferential cross section of p- 3 H scattering. However, the 
resonant behaviour of the nuclear scattering amplitude is 
barely visible in the data. In contrast to the electromag- 
netic probe that amplifies the resonance signal consider- 
ably (see Fig. 1 of Ref. [H). 

Conclusions. We have calculated the isoscalar 
monopole 4 He — > 4 He* transition form factor FmW) 
with realistic nuclear forces (N 3 LO+N 2 LO, AV18+UIX). 
Unexpectedly the results are strongly dependent on 
the Hamiltonian. Therefore this observable is ideal 
for testing nuclear Hamiltonians. As surprising as the 
large potential model dependence, is the fact that our 
Fm results are at variance with the experimental data, 
particularly large differences are found in the case of the 
chiral forces. It is very unlikely that corrections to the 
isoscalar monopole operator can lead to large effects. 
In order to clarify the situation it is highly desirable to 
have a further experimental confirmation of the existing 
data and in particular with increased precision. On the 
theory side further insight could be gained by an analysis 
of sum rules, densities, and different 3NFs. These issues 
will be the subject of future studies. 
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